o 
o 
o 

(N 

Oh: 

D 
00 



> 

in 

m 

o 
o\ 
o 
o 
o 

:a 

Oh: 



X 



Two-photon reactions at high energies 
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We present a short overview over the different contributions to high energy photon photon scattering and 
explore the possibilities of tuning the sizes of the scattered objects by changing the virtuality of the photons. 
We compare the experimental data with model calculations. The difficulties with inclusive charm production are 
discussed. 



1. Introduction 

The reactions 7^*) 7^*) -^ hadrons at high ener- 
gies can be viewed as scattering processes of par- 
tonic systems with tunable size, since the reaction 
time of a highly energetic photon is short as com- 
pared to its dissoziation time into partons [?] and 
the size of the system depends on the virtuality 
of the photon. The situation is sunrmarized in 
figure 1^. 

For high virtualitics of both photons not only 
the parton densities in the photon can be de- 
scribed well by perturbative QCD but also the 
interaction between the small partonic systems. 
Therefore we have in -yy*' j^*> sytems the unique 
possibility to isolate a small angle high energy 
scattering process determined by perturbative 
QCD. This investigation is particularly interest- 
ing since it can clarify the question of the strong 
energy dependence at high photon virtualitics ob- 
served at HERA. In deep inelastic scattering the 
energy dependence at fixed photon virtuality Q^ 
is expressed through the dimensionless Bjorkcn 
variable x = Q^/{Q^ + W^) where W is the 
j{*) ^(*) c ni. energy. There are currently two 
principal explanations for the strong dependence 
on 1/x under discussion. One starts from the 
usual QCD evolution using renormalization group 
arguments to sum powers of a^ \og{Q^/Ql). Here 
the x-dependence cannot be predicted but there 
are initial conditions at fixed Qq which de- 
scribe the observed cross sections very well. The 
rapid rise of the proton structure function with 
small X is then explained by a singularity of the 



Mellin transform of the DGLAP splitting func- 
tion. On the other hand a summation of pow- 
ers of as log(W^^/Wn ) first performed by Lipatov 
and collaborators [g|j|] yields a power behaviour 
of the total cross section like {W'^ /WgY*^ where 
to leading order e/i ~ as41og2 ss 0.55 even for 
a small value a^ = 0.2. If this behaviour is in- 
terpreted in the language of Regge trajectories it 
correponds to a trajectory with a much higher 
intercept than the soft pomeron needed to de- 
scribe high energy hadron hadron scattering Wj. 
The summation of powers of as log(W^^/WQ )is 
a perturbative procedure and the corresponding 
pomeron is called the hard pomeron or from the 
initial of its authors BFKL-pomeron. The BFKL 
formalism has been applied to 7*^*) 7*^*^ scattering 
in |,|. 

The nice picture seemingly explaining the 
strong rise at small x was disturbed by the calcu- 
lations in next-to-leading-order, again by Lipatov 
and Fadin R] and an independent group B. The 
contributions to the intercept where as large as 
the ones from the leading order and reduced the 
intercept to a value very near to one. This clearly 
indicates that resummations are necessary. At- 
tempts in this direction |^,|l^ lead an intercept 
near 1.2. 

Objections have been brought forward against 
the DGLAP evolution, and a solution [|ll| which 
is also very successful from a phenomenological 
point of view has been proposed, namely the in- 
troduction of two pomerons; a soft one relevant 
for low Q^ with intercept 1.08 and a hard one 
important for high Q^ values with intercept 1.42, 
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Figure 1. The diflFractive contribution to hadronic 
^(*) ry{*) scattering. The size of thw quark anti- 
quark system is inversly proportional to the vir- 
tuaUty: Ri ^ l/Q., = y/-qf. 



i.e. near the leading order BFKL intercept. Since 
in 7* 7* reactions with highly virtual photons no 
non-perturbative input as the initial conditions of 
the parton distributions of the proton arc needed 
this reaction plays a priviliged role for the clarifi- 
cation of the nature of the energy dependence in 
perturbative QCD. 

The experiments of LEP2 open thus a new 
window to investigate diffractive prozesses, but 
the kinematic range is not such as to make non- 
diffractive and soft diffractive contibutions negli- 
gible. We put in this contribution the main em- 
phasis on the soft diffractive contributions. They 
are modelled by an approach to high energy scat- 
tering and non-perturbative QCD which turned 
out to be successful in many reactions. The total 
ry{*) ry{*) cross section is in this description given 
as a superposition of the forward scattering am- 
plitude T{W^,Ri, zi, i?2, 2:2) of two colour-singlet 
quark-antiquark pairs |12| with size Ri and R2 
respectively with the p-y the quark densities in a 
photon (see ^ as weight factors 

O'-y'ga* = / (fRl(fR2 / dZldZ2P'y{Ql,Rl,Zl) 

PyiQlR2, Z2)T(M^^ i?i, 21, i?2, 22). (1) 

If in a scattering prozess one object is small and 
one large the variation of the amplitude with the 



size of the small object can be calculated both 
in perturbative QCD and in non-perturbative 
models. They both predict a variation (apart 
from logarithms) of the forward scattering am- 
plitude ~ R^ and hence ^ 1/Q^. Therefore it is 
not astonishing that both perturbative and non- 
perturbative calculations give reasonable results 
also outside their domain of strict applicability. If 
however both objets are of equal size and small, 
i? = i?i « i?2, the behaviour is very different: in 
perturbative QCD the cross section is proporti- 
nal to R^, in the non-perturbative model how- 
ever like i?"*. This opens the aforementioned sep- 
aration of the perturbative domain. But even in 
this case the perturbative model is useful since it 
allows to estimate the kinematical region where 
non-perturbative effects can be neglected. 

As non-diffractive contributions we have the 
box diagram ( figure (a) ) and the reggeon ex- 
change, indicated in figured (b). 
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Figure 2. a: The box diagram giving rise to a 
fixed pole in the complex J plane, b: A schematic 
description of the reggeon exchange. 



2. Non-difFractive and soft diffractive con- 
tributions to 7*^*) 7*^*) scattering 

The box diagram of figure y (a) plays a special 
role. For very large virtualities of the photons 
or high quark masses the QCD corrections to the 
box diagram are under control and it should not 
be modified by them in an essential way. In Regge 
language the box diagram gives rise to a fixed 
double pole at J = in the angular momentum 
plane. Its forward amplitude without any approx- 
imations is given in |l3] . It decreases with energy 



~ \og{W)/W^, it is important at moderate en- 
ergies and large virtualities. 

The reggeon exchange has to be added to the 
flxpole contribution. It is perhaps the biggest 
source of uncertainty in the model calculations. 
We use the factorizing expression |13||: 



7*7* (S, Wli Q2) — 




22'-'/' '^ \1^ 

''''' A^Ql + A^ 


^Ql + A 



(2) 
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with A = 0.3 GeV^, C = 0.38,and 7/ = 0.45. It 
shows the typical energy behavoiour of the ex- 
change of the reggeon trajectory. 

In the limit of one the photon virtualities Q^ 
going to zero it agrees for x > 0.1 very well 
with the hadronic photon structure function re- 
lated to the DGLAP evolved pion structure func- 
tion by naive VDM . It might somewhat over- 
estimate the reggeon contribution at small vir- 
tualities: If it is extrapolated to 77 scattering 
the Regge-contribution to a-y^ is about a factor 3 
larger than expected from factorization and the 
values obtained in |11[ . 

The soft dijfractive contribution is evaluated in 
an approach UM to high energy scattering par- 
ticularly suited to incorporate non-perturbative 
aspects of QCD. The non-perturbative behaviour 
of QCD is treated in the model of the stochas- 
tic vacuum |na,n6[. It is based on the assump- 
tion that the infrared behaviour of QCD can 
be approximated by a Gaussian stochastic pro- 
cess in the gluon field strength. With this as- 
sumption one can evaluate the expectation value 
of Wegner- Wilson loops and one obtains con- 
finement for non-Abelian theories. AUthough 
the model was originally formulated in Euclidean 
space-time it can be extended to describe the ex- 
pectation value of two Wegner- Wilson loops in 
Minkowski space with lightlike sides ||l^,|l^,|l^ . In 
this way one obtains the non-perturbative contri- 
bution to the scattering amplitude of two quark- 
antiquark pairs. 

The cross section for scattering of quark- 
antiquark pairs with sizes i?i and R2 respectively 
obtained from this model can be approximated to 
a reasonable accuracy (10-15 %) by the following 



simple factorizing formula: 

Cdipdip{Rl,Fi2) — 0.7 — —^ — Ri R2X 



(3) 



where {g'^FF) is the gluon condensate [|9[, and a 
the correlation length of the gauge invariant gluon 
field strength correlator, the numerical values are 
given below in equation ^ 

As can be seen from this equation the 
model depends essentially on two typically non- 
perturbative parameters, which specify the Gaus- 
sian process mentioned above: the strength of the 
gluon correlator given through the gluon conden- 
sate and a, its correlation length. From these on 
can calculate the slope of the linear confining po- 
tential [p5 16|. As it stands the model leads to 



cross sections which are constant with increasing 
energy. The parameters of the model were fine 
tuned by a fit to the iso-singlet exchange part of 
(anti-)proton-proton scattering Sit W = \fs~ 20 
GeV. The resulting parameters which are in the 
limits of the values determined by lattice calu- 
lations and also leed to the correct value of the 
slope of the confining potential are: 



{g'^FF) = 2.49 GeV* 



0.35 fm 



(4) 



The phenomenologically observed increase with 
energy of hadronic total cross sections like 
g(Qp-i) T^\^]^ Q,p ~ j^ Qg m can be incorporated in 
two ways: either one lets the radius of the hadrons 
increase with s p7| , p0| -p^ , or one takes the model 
as a determination of the Regge residue and adds 
the Regge-like increase with energy by a factor 
{W^ /WgY"'^-^^ with Wo = 20GeV. These two 
approaches give very similar results, and we adopt 
the latter in this paper as it is the more conve- 
nient in the present context. 

Whereas hadron-hadron scattering and soft 
electroproduction processes (i.e. those with low 
photon virtuality Q^ ) can be described very well 
in this way, this is not the case for hard elec- 
troproduction processes where the energy (1/x) 
dependence is much stronger than indicated by 
the soft non-perturbative pomeron with intercept 
ap — 1.08. As discussed in the Introduction the 



occurence of a second (hard) pomeron as pro- 
posed in |lj] can explain the data in a consistent 
way. This two pomeron approach was adapted 
to the MSV model in ||2^ and very successfully 
tested for the electro- and photoproduction of 
vector mesons and, more relevantly here, for the 
proton structure function over a wide range of x 
and Q^. As in ||ll|] it was found that the soft- 
pomeron contribution to F2, after an initial in- 
crease with increasing Q^, has a broad maximum 
in the region of 5 GeV^ and then decreases as Q^ 
increases further i.e. it exhibits higher-twist like 
behaviour. In the context of the present model 
this is a consequence of the decreasing interac- 
tion strength with decreasing dipole size. 

It is worth recalling the salient features of this 
version of the two-pomeron model which have 
been slightly modified for the theoretical calcula- 
tions presented here, to illustrate the distinction 
between the soft and the hard pomeron in dipole- 
dipole scattering. In |23| it was assumed that all 
dipole amplitudes in which both dipoles are larger 
than the correlation length a = 0.35 fm are dom- 
inated by the soft pomeron, and the energy de- 
pendence therefore given by (W^^/VFq)'^'"-''' with 
Wo = 20 GeV and esoft = 0.08. This ensures 
that the hard pomeron has essentially no impact 
on purely hadronic scattering. If at least one of 
the dipoles is smaller than a — 0.35 fm then the 
energy dependence is replaced by {W^ /WqY'^'^'''' 
with etard = 0.42 . It turned out that the model 
overestimated the non-perturbative contribution 
of very small dipoles so it was put to zero if ei- 
ther of the dipoles is less than 0.16 fm. With 
only four parameters it was possible to obtain a 
good description of data for the proton structure 
function and for the electroproduction of vector 
mesons without noticeably affecting earlier fits to 
hadron-hadron scattering. 

We apply this two-pomeron model to the eval- 
uation of the 7(*)7(*) cross sections. It should 
be noted that the simple factorisation formula 
aj~^ = a'^p/app is no longer applicable in the two- 
pomeron situation. 

The considerations outlined briefly above lead 
to a model for the scattering of quark-diquark 
dipoles on each other approximated by y. In or- 
der to relate it to ^'■-*'>j^*^ interactions we have 



according to |l| to introduce the photon wave func- 
tion. 

The quark-antiquark density in a photon with 
virtuality Q^ = —q^ and helicity A is to lowest 
order perturbation theory given by: 

^e}Q'z^il-zf,Ko{eRr 
p^(A = ±l,Q2,i?,z) = ^^eJ.x 

({z^ + (1 - zf)e^Ki{eRf + m}Ko{eR)A . (5) 

Here if is the charge of the quark in units of the 
elementary charge, i.e ±-i, ±|, nif is the La- 
grangian quark mass; A = indicates a longitudi- 
nal, A = ±1 a transverse photon. The singularity 
of the Bessel functions at i? = does not cause 
any problems, since for the evaluation of observ- 
able amplitudes the density is multiplied by the 
dipole cross section and the latter is proportional 
to R^ at small values of R , see 13 and 0. From 
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Figure 3. The quantity 1000 / aR^ J^ dzp^{X = 
1,Q^ < R,z) as function of R for different values 
of Q^ for light quarks, solid <=> Q^ = 0, long 
dashes <:^ Q^ = 10 GeV2, short dashes 4^ Q"^ = 
100 GeV2, and for charmed quarks, Q-^ = <4> 
dot-dashed line. 



we see that the scale is set by \/z{l — z)Q^ 



For longitudinal photons the factor z^{l — z)"^ in 
the density pj implies that the main contribution 
comes from the region z w i and the relevant 
scale is Q^/4 + m^. For transverse photons how- 
ever the endpoints z — 0, z = \ are not sup- 
pressed and a transverse photon might leak into 
the non-perturbative region even if the virtual- 
ity is quite high. This can be seen directly from 
figure 1^. 

The expressions g are reliable for large values of 
Q^ and/or large quark masses (i.e. c and 6). For 
small values of Q^ and light quarks the distance 
of the quark antiquark pair increases and confine- 
ment effects will become important. Often vector 
meson dominance is applied, another possibility is 
to introduce a constituent quark mass as infrared 
regularisator in the photon wave function. This 
approach will be justified in the following, it al- 
lows a very economical description of the photon 
wave function at low virtualities, which interpo- 
lates smoothly to high virtualities and even takes 
care of the "hard part of the photon" at low vir- 
tualities. This is especially relevant for photon- 
photon interactions. 

Let us begin with a model investigation [g4| 
with scalar photons and quarks. In such a case 
the "photon" -wave function at high virtualities is 
given by: 

~ ^ 1 



t±2 +^(1 _ ;2)Q2 _|_,„2 



^^(R) = _K^(^^z{\-z)Q-^^Tn)R) (6) 

For low Q^ we expect confinement to modify these 
perturbative expressions considerably. 

The structure of || is the same as that of a non- 
relativistic Greens function for the relative mo- 
tion of a free two body system with reduced mass 
m 



Go(i?,0,M) = — i^o(V2TOM|i?|), 



(7) 



where vnM = —mE stands for the virtuality Q^ . 
In order to impose confinement we go from 
the free particle state to a system hold togeter 
by a harmonic oscillator potential. As has been 
pointed out in reference [|9| the harmonic oscil- 
lator is a very useful model for QCD: it shows 



both confinement and asymptotic freedom. We 
therefore investigate the effects of confinement by 
comparing the free Greens function [7| with that of 
the full harmonic oscillator 



Gff(i?,0,A/)= Y. 



{ni + n2 + l)uj + M 



(8) 



which can be calculated easily. 

As has been shown in |^^ a free Greens func- 
tion with a M-dependent shift: M ^ M + s{M) 
yields an excellent fit to the exact Greens func- 
tion ^ for all M > 0. We transfer this procedure 
to QCD by performing a Q^-dependent shift of 
the flavour mass trif —> ruf + m{Q^) in the per- 
turbative photon wave function. The numerical 
value of the shift m{Q'^) can be fixed by a fit 
to the phenomenologically known vector-current 
two-point function and thus no new parameter 
is introduced. The following linear parametriza- 
tions can be used 



muM - 



"^o(l-S) 



Q2 < 1.05 
Q2 > 1.05 



(9) 



mc 



rrir 



0.15 + 0.16(1- ^ 



0.15 



g2 < 1.6 

Q2 > 1.6 



= 1.3 



The parameter ttiq for the u, d quarks was 
found to be mo = 0.21 ± 0.015 GeV.We see 
that above 1.05 for light quarks and 1.6 for 
strange quarks the current mass can be used. For 
charmed quarks the perturbative expression as it 
stands can be taken. 

The densities for transverse photons integrated 
over the longitudinal momentum fraction z and 
multiplied by 1000 R^ are shown in figure g. The 
longitudinal photons are more concentrated at 
small separations R as explained above. 

3. Results 



3.1. Total cross sections | 13 | 

We first discuss the results for the total cross 
section of real 77 scattering. Using the mass 
rriq = 0.21 GeV m, = 0.31 GeV rUc = 1.3 GeV 
we obtain the cross sections given in figure 0. 
The experimental cross section shows an increase 
with energy which is distinctly stronger than for 
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Figure 4. Cross section (nb) for 77 scattering as 
function of the cm. energy W (GeV) compared 
with L3 data. L3: Triangles |2^] The soUd curve 
is our model without adjusted parameters. It con- 
sists of the following contributions: soft pomeron 
<^ long dashes; hard pomeron <^ short dashes; 
fixed pole (box) <^ dot-dashes; reggeon <(4> dots. 



hadron hadron scattering. Our model reproduces 
this result and it can be explained very natu- 
rally. The singularity of the photon density p^ 
at i? = which reflects the pointlike coupling 
of the thoton to the quark-antiquark pair gives in 
7 7 scattering a stronger weight to smaller dipoles 
in hadron hadron scattering. Since small dipoles 
couple to the hard pomeron which has a consider- 
ably larger intercept we obtain an overall energy 
dependence which is stronger than for hadrons. 
The contribution of the box diagram (dot-dashes 
in figure 0) and the reggeon exchange (dots) are 
of minor importance above 14^ = 30 GeV. As ex- 
pected from the discussion in the introduction the 
purely non-perturbative model decreases much to 
fast with virtuality if both photons have approx- 
imatively equal values for Q^ . This can be seen 
in figure ^ where the Q^ dependence of the to- 



tal 7^*) 7*^' cross section is diplayed as a func- 
tion of Q^ = {Q\) ~ (Ql) at the fixed cm. en- 
ergy W Ki 20 GeV; at Q^ = 14 GeV^ the non- 
perturbative contribution is negligible. The en- 
ergy dependence of the 7^*^ 7^*^ cross section at 
Q^ = 3.5 and 14 GeV^ is shown in figure ^. If the 



Figure 5. Cross section (nb) for 7*7* scat- 
tering as function of the average photon virtu- 
ality (Q2) ^ {qD ^ {qD (GeV^) a cm. energy 



ly « 20 GeV compared with L3 data, 
as in figure 
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theoretical non-perturbative and non-diffractive 
contributions are subtracted from the experimen- 
tal values the resulting energy dependence can be 
very well described by a power behaviour (W^Y 
with 6 = 0.3... 0.4. 



3.2. Photon structure functions [13| 

In a photon strucurc function one photon is real 
and the other one has a virtuality Q^ ^ 0. The 



structure function F2 is defined as 



F^{x,Q') 



Q' 



'{W,Q^) 



(10) 



with X — Q^/{W^ + Q^). As mentioned in 
the introduction in that case both perturbative 
and non-perturbative QCD give the same Q^ be- 
haviour and we see indeed in figure H that the 
non-perturbative model describes the data very 
well up to Q^ — 20 GeV^. In order to extrapolate 
nontrivially to Q^ = a modified structure func- 
tion F^ = 9l±o^F^(^x = ^^,Q^) has 
been displayed and compared with the correpond- 
ingly multiplied experimental data from different 
experiments. In figure ^ the photon structure 
function for Q^ — 10 GeV^ is shown. The agree- 
ment between theory and experiment is similarly 
good for other virtualities. 
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Figure 6. Cross section (nb) for 77 scattering 
as function of the cm. energy W (GeV) com- 
pared with L3 data for (Q^) = 3.5 and 14 GeV^ 
respectively. Symbols as in figure 
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Figure 7. The modified photon structure func- 
tion Fl = (Q2 + o.6GeV2)/Q2i^2'(a; = (3V(W^^ + 
Q^),Q^) as function of Q^ for W ^ 20 GeV. 
The data ^^^ are: Triangles -^ L3 ; Dia- 
monds,Boxes <^ OPAL ; Stars ^ ALEPH The 
solid curve is our fit without adjusted parame- 
ters. For the different contributions see caption 
of figure |[ 
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Figure 8. The photon structure function 
F'^ix, Q"^) as function of x = w'^+o'^ ^°'' Q^ — ^^ 
GeV^ . The data H-Q are: Triangles ^ L3 ; 
Diamonds 4^ OPAL ; Stars ^ ALEPH The solid 
curve is our fit without adjusted parameters. For 
the different contributions see caption of figure 0. 



3.3. The reaction 7**' 7 -> ccX 
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Figure 9. The cross section for the reaction 
77 ^ ccX . The sohd Hne is the direct contri- 
bution to the structure function as predicted by 
the model from a picture correponding to figure 
and the box digram with charmed internal quark 
lines (dot-dashed hne); data points are from |28| . 
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Figure 10. The photon structure function _FIJ^ at 
(Q2) = 20 GeV^. The solid hne is the direct con- 
tribution to the structure function as predicted by 
the model from a picture correponding to figure |l| 
and the box digram with charmed internal quark 
lines (dot-dashed line). The boxes are the data 
for the charm structure function |g^, the trian- 
gles are the results for the full structure function 
at the same average Q^. 



New measurements from L3 ||2§] and OPAL Q 
yield data for inclusive charm production. In 
figure 1^ the total cross section for the reaction 
77 ^- ccX is given. We have calculated diffrac- 
tive charm production according to the physical 
picture of figure |^ where the quarks in the loop 
of photon 1 are charmed, in loop 2 all flavours 
are summed over and correspondingly with 1 
and 2 interchanged; for the box diagram only 
charmed quarks contribute. We see that the the- 
oretical values for the direct diffractive and non- 
diffractive production underestimate the data ||2^ 
by more than a factor two. The same happens for 
the case if one photon is virtual. In figure |l^ the 
charm contribution to the photon structure func- 
tion is displayed for (Q^) = 20 GeV^. The central 
value of the experimental point ||2^ at (x) = 0.02 
is again far above the theoretical value. 

In principle this could indicate a failure of the 
model prediction for the cross section || for the 
case that one object is small and one large, but 
this is unlikely for the following reasons. 



1) The charm structure function of the proton IQ 
agrees very well with the predictions of the model 
as can be seen from figure 11 and here we have 



the same situation of one object large and one 
small. 

2) The full photon structure function at Q^ = 10 
GeV^ is correctly predicted by the model, as can 
be seen from figure |[ At Q'^ — 10 GeV^ the 
light quark density in a photon peaks at about 
the same quark- ant iquark distance as a charmed 
pair in a real photon (see figure ^. 

So we would like to conclude that there must 
be other mechanisms of charm production but re- 
frain from speculations until the experimental sit- 
uation especially of the charmed photon structure 
function is clearer. 
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Figure 11. Cross section for the reaction 7*p ^ 
ccX. The solid curve is our model for the di- 
rect contribution for photons with Q^ = 10 GeV^. 
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Discussion 

S. Narison, LPMT: / am, a hit ennoyed with 
the log(M^/rn^) behaviour of the cross section 
when ruq ^0. / would have expected a be- 
haviour like \og(W/ ^) where ^ is a subtraction 
point. What is the meaning of rriq for light quarks 
Q 

H.G. Dosch: For real photons the box diagram 
has an infrared singularity which is regularized 
by the quark mass. The same happens for the 
pomeron exchange of real photon photon scatter- 
ing. Here the IR singularity leads to 1/m^ ! For 
light quarks an IR regularization procedure has 
to be introduced which at the present state of 
the art is model dependent. 

A. Brandenburg: Could you explain in more 
detail how you obtained the quark mass that regu- 
larized the photon wave function? What is the 
physical interpretation? Does it have anything 
to do with the mass parameter in the QCD La- 
grangian ? 

H.G. Dosch: The regularization of the photon 
wave function by an effective mass was motivated 



by model investigations [gj]. The values for the 
regularizing masses for light quarks are given in 
equation 9 and are close to " constituent masses" . 
They were obtained in the following way: The 
second derivative of the lowest order two point 
function for the vector current was compared with 
the corresponding phenomonological expression, 
consisting of the p— pole and a continuum. The 
theoretical expression contains the quark mass as 
only free parameter and it was fixed by equating 
the two expressions. For Q^ > 1 the agreement 
could not be improved by a finite quark mass and 
there the Lagrangian mass can be taken, but it 
is of little influence anyhow. For heavy quarks 
the Lagrangian quark mass at the scale Q^/4 is a 
reliable regularizator. For a detailed discussion I 
refer to reference EJ]. 
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